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Abstract 

We study Delbriick scattering in a Coulomb field in the presence of a laser field. The amplitudes 
are calculated in the Born approximation with respect to the Coulomb field and exactly in the 
parameters of the laser field having arbitrary strength, spectral content and polarization. The 
case of high energy initial photon energy is investigated in detail for a monochromatic circularly 
polarized laser field. It is shown that the angular distribution of the process substantially differs 
from that for Delbriick scattering in a pure Coulomb field. The value of the cross section under 
discussion may exceed the latter at realistic laser parameters that essentially simplify the possibility 
of the experimental observation of the phenomenon. The effect of high order terms in the quantum 
intensity parameter x of the laser field is found to be very important already at relatively small X- 
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I. INTRODUCTION 



The theoretical and experimental investigation of nonlinear Quantum Electrodynamics 
(QED) effects arising due to vacuum polarization in the presence of a strong external field 
renders possible testing QED when the nonperturbative effects in the field strength (con- 
tributions of higher-order terms) are very significant. One of these effects, the scattering 
of a photon in the electric field of atoms (Delbriick scattering), has been investigated in 
detail both theoretically and experimentally . Also, a decisive progress has been 

obtained in the understanding of the process of photon splitting in an atomic electric field 
0, 0]. In [?| the first successful experimental observation of photon splitting in an atomic 
field is reported. As a result, it has been observed that higher orders of the perturbation 
theory with respect to the parameter Za play an important role and essentially modify the 
cross sections of Delbriick scattering and photon splitting in an atomic field at high photon 
energy (here, Z is the nuclear charge number, a = e 2 = 1/137 is the fine-structure constant, 
with e being the electron charge and the system of units h = c = 1 is used). 

Besides the atomic field, photon splitting has also been investigated in other configura- 
tions of the external field, namely in constant electromagnetic fields (si, [h], 11. 12| and 



in a laser field 131. In 



131 ] photon splitting in a laser field with arbitrary strength and fre- 



quency content has been investigated for an arbitrary energy of the incoming photon. The 
photon splitting amplitude in a laser field depends on two parameters that, in the frame 
where the laser wave propagates in the direction anti-parallel to the momentum k\ of the 
initial photon, are given by 



U x Uq LO X E 

V= a"' X=— TT- (1) 



Here, u)\ is the energy of the incoming photon, uq and E are the characteristic frequency 
and electric field strength of the laser wave, E c = m 2 /\e\ = 1.3 x 10 16 V/cm is the so-called 
critical electric field and m is the electron mass. At present, the parameters of optical laser 
fields as well as of photon sources permit achieving nowadays values of x ~ 1 where non- 
perturbative effects in the laser field become important [ijj]. We emphasize, however, that 
the structure of the photon splitting cross section in the case of both constant and laser 
fields does not allow one measuring this effect without considerable effort. 

In the present paper, we consider a process which is closely related to photon splitting in a 
laser field. Namely, we study photon scattering (Delbriick scattering) in combined Coulomb 



2 



and laser fields. We perform our calculation in the Born approximation with respect to the 
Coulomb field and exactly in the parameters of the laser field. Thus, the amplitude obtained 
is linear in Zot and corresponds to photon splitting in a laser field when an initial real photon 
with the momentum k± = (ui, k±), k\ = 0, splits into a real photon with the momentum 
k% = (k>2, k 2 ), k\ = 0, and a virtual photon (Coulomb quantum) with the momentum 
&3 = (0, q), k 2 = — q 2 . We consider the most interesting case where a nucleus (source of the 
Coulomb field) is at rest and a plane wave, corresponding to the laser field, propagates in the 
direction x anti-parallel to the momentum h\. The plane wave is described by an arbitrary 
four-vector potential A^{ip) = (0, A(0)), where <fi = xx, with x M = (l,x), x 2 = 0, and 
x ■ A(d>) = 0. The amplitudes are derived following the method described in detail in our 



paper 



13j and based on the operator technique developed in |l5|, [l6j for the calculation of 
the polarization and mass operators in a laser field. The advantage of the operator technique 
is that it enables one to calculate the amplitudes of the processes without employing the 
explicit form of the electron propagator in the external field. We first derive the explicit 
expressions of the amplitudes for a general vector potential of the form A(<f)) (Section II). 
Then, we consider the case of a monochromatic laser field with circular polarization (Section 
III) and analyze in detail the properties of the amplitudes for the most interesting parameter 
region: u <^ m and uj\ ^> m (Section IV). Finally, in Section V, the results are discussed 
and the main conclusions of the paper are presented. 

II. GENERAL FORM OF THE AMPLITUDES 

The amplitude M of Delbriick scattering in combined Coulomb and laser fields calculated 
in the Born approximation with respect to the Coulomb field is represented by the Feynman 
diagram shown in Fig. [TJ It has the form 

M = - ^ 2 Ze] - M^e*^, (2) 

where and are the initial and final photon polarization four-vectors, respectively, 
e 3/x = fiVo is the polarization four- vector of the Coulomb quantum, and the formal expression 
for the tensor M.^ yp in the Furry representation reads 

M^ p = ie z [ d A xTr(x\e- lklX Y-~ e lk2X Y e ik:iX Y~~ \x) 

J P — m + iO P — m + iO P — m + iO 

+(k% k£ , v p) , (3) 



Here P = -P /t 7 M , with P^ = id^ — eA^(<f)) and 7^ being the Dirac matrices. 




k 3 

FIG. 1: Feynman diagram for the Delbriick scattering amplitude in combined Coulomb and laser 
fields calculated in the Born approximation with respect to the Coulomb field. This diagram cor- 
responds to the Furry representation and the thick line denotes the electron propagator (Green's 
function) in the laser field. The wavy lines symbolize the external photons, k± and fc^ are the mo- 
menta of the initial and final photons, respectively, fcg is the momentum of the Coulomb quantum. 
The diagram with the permutation kr, «-> fcg has to be added. 



By using gauge invariance it is possible to simplify the calculation of the amplitude M. 
To that purpose we introduce for each photon with the four-momentum fcf (i — 1, 2, 3) two 
four- vectors s^ x ., (Aj = 1,2) 

a° x = 0, (a A ) 2 = l, oi-o 2 = 0, o A -x=0. (4) 
For the virtual photon (Coulomb quantum), we also introduce the third four- vector 

which is perpendicular to k$, e%i, and £32 (see Eq. Then, gauge invariance allows us 

to write the tensor M? vp introduced in Eq. (j2j) in the form: 

M^ = - R ^AxAAm- ( 6 ) 

Ai,A 2 =l,2 
A 3 =l,2,3 
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From this equation we have 

By employing the coefficients i?A 1 A 2 A 3 , we write the amplitude M as 

M = W li Z e) ^ 1 A 2 A 3 (£ 1 ,A 1 e 1 )(e 2) A 2 e^)(£3 J A3^). (8) 

9 Ai,A 2 =l,2 
A 3 =l,2,3 

In this equation it is not necessary to perform any subtractions since it automatically fulfills 
gauge invariance. 

It is convenient to perform the calculations in terms of the helicity amplitudes M CTl0 - 2 , 
with (jj = ±1. In this case, the polarization vector e a of each external real photon with the 
wave vector k and the energy u> satisfies the relations e a ■ k = and e a xk = iauje a . Due to 
momentum conservation in the laser field, the final transverse momenta k 2 ± and k 3 ± (with 
respect to the laser propagation direction) are equal in modulus and opposite in direction: 
k 2 ± = —k 3 ± and k 2 ± ■ x = k 2 ± • ki = 0. We direct the vector a\ along the vector k 2 ± and 
set a 2 to be a 2 = x x a±, then 

ei j<71 = --j={a 2 + latax) , 

e 2 ,<7 2 = -^={a 2 - ia 2 a 2 x k 2 /u 2 ) . (9) 



four- vectors 



The polarization four- vectors e^ CT . = (0,e i>cr J and have the following products with the 

A 4 . 

£l 2,\2 e *2,a 2 = -^(*0-2^A 2) l - 5a 2 , 2 ) , 



esMel = 4a 3 = <W + ^3,3 • (10) 

Using these relations we obtain from Eq. (jSJ) 

M ai a 2 = ( f e) (cria 2 Kn + ^ 22 + io{R l2 - io 2 K 2l ) , (11) 
Zq 



where 



= 1 + " V (^F^ 2 3 • (12) 



(2) 

The calculation of both terms, R\ \ 2 \ an d R 

formed in [3] and we present here only the final results. For the term R^\ 2 \ 3 with A3 = 1, 2, 
we obtain 



? (3) 



, ± l Ai A2A3 a±1KX -"'A1A2A3' * s Q u ite similar to the calculation per- 
ls) 



R 



(2) 

A1A2A3 



me 



(27C) 



■ £ « 

j=l,2,3 



(J) 



(si + s 2 ) : 



■ exp 



-Z'(S1 + S 2 ) - Z — Q 2 ^j,3 



1 exp ^ zS20 + z - 2 y?j _ H s i + s 2j 



S\ + s 2 
1 







x 



2a A , • / dyA( Wj ) + ( Sl + S ^ a Xj ■ Afa) 
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In this expression we introduced the quantities 



(13) 



^(1) — ^A 2 ,A 3 , °(2) — °\i,\ 3 ) (3) — 0Ai,A 2 ) v \ 



5m = 5\ 



"2,3 



2.3 



Q 2 = ^, 

— 1 z 



ifj = l —L- p. ; A(w) = — + 4wu;i/m 2 ) — A{(j))\ , Q = u>2 — 0J\ ■ 

S\ + S 2 Tfl 



(14) 



(2) 

Note that R\^x 2 x 3 contains only odd powers of the external field. 

We pass now to the term R^_^^ 2 ^. In order to write it in a more compact form, we 
introduce the convenient notation 



c . . s 2 v 2 - S3V3 

= A(s 2 r 2 y + s 3 r 3 (l - y)) , D 2 (y) = A(s 2 r 2 y) , D 3 (y) = A(s 3 r 3 y) 
1 Z" 1 

F =Sj d yisiDi(y) + s 2 D 2 (y) + s 3 D 3 (y)}, 

F = SF 2 - I dy[ Sl D\{y) + s 2 D 2 2 (y) + s 3 D 2 3 (y)] , 



G = ^ /" dy [saCACy) - Dad/)) + ^(y) - D 3 (y))] 



S 



k? 



m 



JVi = F-^, iV 2 = F-D 2 (l) + ^ 

5 ' v ; Sz/ 2 



iV 3 = F - D 3 (l) + ^ , V = i(di + cZ 2 + d 3 ) , 



d 2 (i)-d 3 (i; 

2ri 



2r 9 



2T3 



(15) 
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We obtain for A 3 = 1,2: 

oo 

(3) ime 3 ^ fff dsids 2 ds 3 



R 



AiA 2 A 3 2yj-2 



V 



s 2 



d(f> 



x exp 



^i^s 2 s 3 s lS3 2 

+F -S + P -G 



x ^ - 4(a Al • iVi)(a Aa • iV 2 )(a A3 • iV 3 ) + 4(a A2 • iV 2 )(a A3 • N 3 )(a Xl ■ V) 
-4v 3 (a Xl ■ NJia^ ■ N 2 )(a x . A ■ V) - 4z/ 2 (a Al • A^K, ■ iV 3 )(a A2 ■ V) 
+(a Al • JVi)! ^^ + + [p • (di - da - d 3 )] A2A3 - 2i Wl [d 1 • (d 2 + d 3 )] A2As 



-{ax* ■ N 2 ) 



1v 2 v 3 2v 3 

[P P]AlA3 + ^(2) + [P ' (di + d 2 - d 3 )] AlAs - 2^ 3 [d 3 • (d x + d 2 )] AlAa 
lv 3 lv 3 



-(a A3 • Ar 3 )| [P 2 ^ AlA2 + [p • - d 2 + d 3 )] AlA2 - 2z/ 2 [d 2 • (di + d 3 )] AlA2 
+(a Al • p)[p • di] A2As - — (a A2 • p)[p • d 3 ] AlA3 - — (a As • p)[p • d 2 ] AlA2 

^2 ^3 

j [ a Ai • di] A2A3 [a A2 • d 3 ] AlAs [a Ag ■ d 2 ] AlA2 1 |^-[a A3 • d 2 ] AlA2 

A I V 2 V 3 ) ZV 3 



+ 



^2^3 

Q 2 



-(a Al -p)(a A2 -p)(a A3 ■ p) 



4z/ 2 z/ 3 



[5 ( i) (a Al • p) + 5 (2) (a A2 • p) + (a As • p)\ 



+ ^7^(l)(«A 1 ■ P) + 5(2) (flA 2 • P) - 5(3) («A 3 ' P)] 

2i 

+ - [5 ( i } a Al • (ATi - V) + 5 (2) a A2 • (iV 2 + v 2 V) + 5 (3) a A3 • (N 3 + v 3 V)\ 



+Ui -v 2 XJ 2 - v 3 U 3 + 2u 2 u 3 U 4 



V 



(16) 



The operator V in front of the integral takes the odd part of the integral with respect to 
the external vector potential We have also introduced the operator 

[X ■ Y] XiXj = (X ■ a Xi )(Y ■ a X] ) + (X ■ a Xj )(Y ■ a Xi ) - (X ■ Y)(a Xt ■ a X] ) (17) 
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for any two vectors X and Y and we employed the abbreviations 

U x = -(a Xl ■ p)[d 2 ■ d 3 ]x 2 \ 3 + (a Aa ■ p)[d 2 ■ d 3 ] Xl \ 3 + (a Ag ■ p)[d 2 ■ d 3 ] Xl \ 2 
+2(a Al • p)(a x , 2 ■ d 2 ){a X3 ■ d 3 ) - (d 2 ■ p)[d 3 ■ a A2 ] AlAs - (d 3 ■ p)[d 2 ■ a As ] AlA2 , 

U 2 = (a Xl ■ p)[dx ■ d 2 ]A 2 A 3 - (a X2 ■ p)[d x ■ d 2 ] AlAs + (a Ag • p)[d\ ■ d 2 ] XlX , 2 
+2(«a 2 • p)(oai • d 2 ){a X3 ■ d x ) - (d x • p)[d 2 ■ a X3 } XlX2 - (d 2 ■ p)[d x ■ a Al ] AaA3 , 

U 3 = (a Xl ■ p)[di ■ d 3 ]A 2 A 3 + («a 2 ■ p)[di ■ d 3 ] AlAs - (a Ag ■ p)[d x ■ d 3 ] XlX2 
+2(a As • p)(a Xl ■ d 3 )(a X2 ■ d x ) - (d x ■ p)[d 3 ■ a X2 } XlXs - (d 3 • p)[d x ■ a Al ] A2Aa , 

U 4 = (a Al • d 2 )[di ■ d 3 ] A2Aa + (a Al ■ d 3 )[di ■ d 2 ] X2X[i - (<*Ai ■ d x )[d 2 ■ d 3 ] AaA3 

+2(a Al • d x )(a x , 2 ■ d 2 )(a x , A ■ d 3 ) + (d x ■ d 2 )(d 3 ■ {5^a X2 - 5( 3) a x J) 

-(d x ■ d 3 )(d 2 ■ (5 {2) a X2 - 5 (3) a A3 )) - (d 2 ■ d 3 )(di • (5 {2) a X2 + 5 {3} a X3 )) . (18) 

For A 3 = 3 we have 



x U(t x + r 3 )(a Al • iV!)(a A2 ■ N 2 ) 

+2(a Al • ATi) [4y 2 n(a A2 • V) - p ■ a X2 + 2z/ 2 z/ 3 (a A2 ■ (d 2 + d 3 ))] 
-2(a A2 • iV 2 )[4ri(a Al • V) - p • a Al + 2z/ 3 (a Al ■ d 3 )] 
-4z/ 2 r 3 [d 2 • d 3 ] AlA2 - 4^ 2 Ti[d! ■ d 2 ] 

AiA 2 

+2rx[p • di] AlA2 + 2r 2 [p • d 2 ] AlAa + 2r 3 [p • d 3 ] XlX , 2 



Equations ( |T3l) . ( |T6l) and ( |T9l) are valid for an arbitrary vector potential A(</>). The structure 
of Eqs. f|T3|) . (fT6l) and f|T9|) is very similar to those obtained in the investigation of photon 

8 






(19) 



splitting in a laser field. The main difference is the presence of the term proportional to Q 2 
in the exponent. However, as it will be seen below, this term is crucially important for the 
properties of the amplitudes. 



III. THE CASE OF A MONOCHROMATIC, CIRCULARLY POLARIZED PLANE 
WAVE 

We consider now the laser field to be a monochromatic plane wave with frequency 
Then, after the integration over <fr in Eqs. f fT3|) . f|T6|) and f|T9|) . we arrive at the following 
form for the helicity amplitude M aia2 : 



M, 



<T 1 (T2 



y j 2ir5{u 2 nu )M njai(T2 . (20) 



Here the prime indicates that the summation is performed over odd numbers n, and M nj(TlCT2 
is the invariant amplitude of the process with absorption of n laser photons if n is positive or 
emission of |n| laser photons to the laser field if n is negative. Unlike the process of photon 
splitting into two real photons, in our case n can be negative with the only restriction 
n > —uji/ujo, following from the condition uo 2 > 0. 

The cross section corresponding to the amplitude M njCT1(J2 has the usual form 



2 co»2 sin 9 d9 dip 



da n a , a2 \M n axC! 2 1 „ 1 (21) 

where 9 is the angle between vectors &2 and ki, and if is the azimuth angle of the vector 
k 2 ± in the plane perpendicular to X. In terms of Q = uo 2 — 0J\ and 9 we have 

xk 2 = {u)\ + fi) (1 + cos 9) , xfc 3 = (uji — Q) — (oji + Q) cos 9 , 

\k 2± \ = {ut + tt)sm9, kl = -q 2 = -2^(1 - cos fl) - 2fi 2 (l + cos^) . (22) 

The vector potential of a monochromatic, elliptically polarized laser field is 

A(4>) = A 1 cos(cj o 0) + A 2 sin(cu o 0) , (23) 

with Ai ■ A 2 = Ai ■ x = A 2 • x = 0. In general, the vectors Ai and A 2 are not parallel 
to the unit vectors and a 2 in Eq. respectively, and there is an azimuth asymmetry 
in the cross section for an elliptically polarized laser field (we call to mind that a\ has been 
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chosen to be parallel to k 2 ±). However, in a circularly polarized laser field which we consider 
below, the azimuth asymmetry is absent and the vector A\ can be directed along a\ and 
A 2 along a 2 (for positive helicity (p.h.)) or A 2 along — a 2 (for negative helicity (n. h.)). 
Here, we consider only the case of positive helicity of the laser field while the amplitudes for 
negative helicity can be obtained by means of the relation M aiU2 (n. h.) = M 5 . l52 (p. h.), with 
&i denoting the helicity opposite to <7j. Also, we set \Ai\ = \A 2 \ = A and £ = eA/m. Using 



Eqs. (jTTj) and (1121 . we represent M nj 



as 



n,<Ti<T2 



Af (2) + M (3) 



(24) 



where for M^l ia2 we find from Eq. ffl~3 



n,H — 

r(2) 
'n, — h 

r(2) 



-^(G^ + G^); 
5n,-i(G[ +) + G 3 +) ) + + Gff >) ; 



M (2 | 
M, (2) 

u 3 uiQ 2 J s J 

sin t? \ sin $ 
+ 



xe 



id 



2i I e i,? - 



1 + w(l - u)Q% t3 + £ 2 ( 1 - 

1? = ±2tt(l — u)S7]Vj . 



sin 

tf 2 



•d J w(l — w 

(3) 

In order to write Mn,a 1 a 2 in a compact form, we introduce the following definitions: 

e 



(25) 



c 



TlV \S 2 T 2 S 3 T 3 



>) = arg C 



* = eici 



4rfSriT 2 T 3 \ n 

P S 2 S 3 



^ 3 iAi 2 + % 2 | 2 -% 3 r)-e 2 s ; 



r 2 



T 3 



5 



5 = Si + s 2 + s 3 . 



In terms of these functions, the contribution MA,^^ has the form 



(26) 



iZa 2 p f T T dsi ds 2 ds 3 

^hQiJJJ S 2 



3 

J=-3 



exp 



i Un V S -^ - ^?g 2 + F - S 



s 



s 



(27) 



where Ji{z) are ordinary Bessel functions. The coefficients Bj >rTl(T2 are expressed via the 
functions fj introduced in Eq. fl26|) and they are presented in the Appendix. The expressions 
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(|25|) and (J27J) are valid for any value of the parameters 77 and £. Note that the expansion of 
the amplitude M n ^ ia2 for a fixed n, contains in general all odd powers k of the parameter 
£, starting from k = \n\. The terms with k > \n\ in M no - lcr2 correspond to rescattering 
processes of absorption and emission of laser photons with a net emission or absorption of 
|n| laser photons. 



IV. ASYMPTOTICS OF THE AMPLITUDES 



In this Section we consider some asymptotics of the amplitudes M„^ lff2 in Eq. 
For optical laser fields we have ujq *C m. Also, at io\ < m the cross section is too small 
in comparison with that of usual Delbriick scattering in a pure Coulomb field. Therefore, 
we consider below the most interesting case, from an experimental point of view, uj\ ^> m. 
In this case, the main contribution to the total cross section comes from scattering angles 



6 < ra/uji <C 1. We also assume that 9 is not too small, namely 6 ^> 6 = ^Jujq/uji. In these 
approximations we have 

^ = %, P = — , Q 2 = p 2 (l + u 3 ). (28) 
4 m 

Note that the amplitudes in the region 9 ~ #0 can also be easily derived from the general 
expressions ( |25|) and (1271) . However, it is difficult to investigate experimentally this region 
of very small 6 and then we do not consider it here. 

We start with the case of small rj = uoquoi/itl 2 and large £ but fixed x = v£ an d P- Nowa- 
days, these conditions can be provided by using available photon sources and strong optical 
lasers [14J]. Note that this limit of parameters provides the fulfillment of the quasiclassical 
approximation. It is convenient to present the cross section (1211) in the form 



dcr n>ai<T2 = cr \ A n ^ 1(J2 \ 2 <j = -a 2 {Za) 2 \ 2 c x 2 , (29) 

P 7T 



where Ac = 1/m is the Compton wavelength and the dimensionless amplitudes A n>++ and 
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A„ i+ _ read: 



l n, ++ 



(S n ,i + / ds / due'^ [1 - 2u(l - u)] 

P Jo Jo 

2 



poc pi pi 2 

+i ds duu dye- l ^y2(-l) n+] J n+3 (z)B j:++1 

Jo Jo Jo j = _ 2 

poo pi pi 

i ds duu dye' 14 ' 1 ^ (-l) n+J ' J n+j {z) B j>+ - . 

Jo Jo Jo „•_ n 



(30) 



Here, the functions ^o,i, -2 an d the coefficients £>j iCri(72 are 

= s + ^m 2 (1 -w)Vs 3 , ^i =^o + 2/(l ~y)u 2 p 2 s, z = Axps 2 u 2 {l -u)y{l-y) , 



£ 2 , ++ = 2 X sV(l - u)[l - 2u(l - u)}(l - 2y) 2 , 



Bi >++ =[l-2ii(l-u)] 



— - + Aipsu 2 y(l - y) 



— ipsu(2 — 3u + 2u 2 ) , 



1 — u fp 1 

B 0<++ = 2i h £-u[-l + 4u(l - u)y(l - y)] - 8xs 2 u(l - u) 2 (l - u + u 2 ) , 



B_ li++ = [l-2u(l-u)] 



h 4ipsu 2 y(l - y) 

P 



+ ipsw(2 - 5m + 6m 2 - 4m 3 ) , 



2p 2 



Boh 



2,++ , 

u 2 (l-u)y(l-y), 



= 4ipsu 2 (l - u) 2 , 
B- 2 ,+- = 8xsV(l- M ) 3 . 



(31) 



The amplitudes A nj and A ni _ + are obtained from A Ht++ and by the substitution 



"71,- 



-A_ r 



(32) 



For uji ^> m it is difficult to measure the energy U2 of the final photon with an accuracy 
of the order of u (energy of optical lasers photons). Besides, it is as well very difficult to 
measure the polarization of the final photons. Therefore, we consider below the cross section 
da summed up over n and over the polarization of the final photon and averaged over the 
polarization of the initial photon, i. e. 



da = - ^2 da nr 



(33) 



n,cri,<72 
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If p <C 1 and x < 1 than the leading terms of the amplitudes are linear in p and the 
nonzero A n ^ ai(72 are 

A n>++ = p(6 n> i - S n ,-1)^1 , = p{5 n ,i^2 + ^,3-^3) , 

2 /*oo /»1 

Fi = - ds due-^°su 2 (l-u)[-3 + Au-Au 2 + Aix 2 s 3 u 2 (l-u) 2 (l-u + u 2 )], 

3 Jo Jo 

4 r 00 r 1 

T 2 = -\ ds due-^°su 3 {l-u) 2 [-3 + 2ix 2 s 3 u 2 {l-uf 
3 Jo Jo 



\2l 

) 



o poo pi 

jr 3 = _,° 2/ ds / d ue -^>sV(l -m) 4 . (34) 
3 Jo io 

Therefore, at p <C 1 

rfa = a (2|J- 1 | 2 +|^ 2 | 2 + |J- 3 | 2 )^, (35) 

P 

and the integrated cross section contains a logarithmic divergence at small p. This divergence 
is eliminated by including the atomic form factor into Eq. ( 1291) and by replacing p 4 in Eq. 
(1291) by Q 4 ~ (p 2 + 4f2 2 /m 2 ) 2 with f2 = nujQ. Thus, we obtain with the logarithmic accuracy 

a = (T (2|J zr i| 2 + l^l 2 + l-T^I 2 ) log ( -3— J , Pmm = max { — !— , — 1 , (36) 

\PminJ [mr scr m) 

where r scr ~ Z~ x ' 3 j{ma) is the screening radius corresponding to the atomic form factor. 
Note that it is also not difficult to obtain the integrated cross section beyond the logarithmic 
accuracy. In Eq. (|36|) the functions Ti depend only on x with J-\(x = 0) = 11/90, T^ix = 
0) = 1/15, and T^{x = 0) = 0. This dependence reflects the contributions of high order 
terms with respect to the laser field. 

We consider now the asymptotics of the amplitudes in Eq. (1301) at p ^> 1 and x > 1. The 
leading contribution to the cross section da comes from n = ±1 so that 

1 1 

A n ,++ = ~\^ n > 1 ~~ bn-l)F 1 , = ~8n,lJ~& j 

poo pi 

Ta = % \ ds du e~ #0 (1 - 2u + 2u 2 ) , 
Jo Jo 

poo pi 

T h = Ai j ds j due'^uil-u). (37) 

Therefore, at p ^> 1 



Jo 



do = oo{2\F A \ 2 + \F b \ 2 ) d ^. (38) 
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Finally, we also present the amplitudes at x 1 and fixed p. This asymptotics corresponds 
to the Born approximation with respect to the laser field. The amplitudes linear in \ are 



Ai,, 
A-i,, 



P\ 'dy f d ^{2y{l-y)u{l-2u + 2u 2 ) 



o 



n 



V 



+— [4y(l - y)u 3 + 2 - 5m + 6m 2 - 4m 3 ] 
-^j2m 2 (1 - u)y{l - y)p 2 [l - 4«(1 - u)y(l - y)\ 

P f dy I ^-{2y{l-y)u{l-2u + 2u 2 ' 



+ - vM 2 - 4m + 3m 2 ) - 2 + 3m - 2m 2 ] 

+^2m 2 (1 - u)y(l - y)p 2 [l - 4m(1 - u)y{\ - y)] 



f 1 , ^duuHl-u) 2 
4p / dy 

Jo Jo 
»1 

.3 



V 2 



2p 2 u 2 y 2 (l- y y 
V 



ip~ I dy 
'o 



1 rfMM 5 (l -m)V(1 -yf 



P 3 



V = l + y(l-y)u 2 p 2 . 



(39) 



Note that the above amplitudes for n — — 1 are, up to an overall factor, nothing but the 
amplitudes of photon splitting in a Coulomb field in the Born approximation. We have 



checked that our amplitudes are in agreement with the results obtained in 



6]. 



V. DISCUSSION OF THE RESULTS AND CONCLUSION 

At iO\ ^> uq, the main background to the process under consideration is the usual Delbriick 
scattering (we denote the cross section of this process as dan) which takes place without 
any laser field. In this case u 2 — oji, and the amplitude corresponds to the exchange of an 
even number of Coulomb quanta. In our case, instead, uj 2 ^ uj\ and therefore there is no 
interference with the usual Delbriick amplitude. It is important that the dependence of da 
and da D on p is very different at small as well as at large p, see Refs. j^, Q]. It is natural 
to compare the cross sections da and da^ obtained both in the Born approximation. In our 
case, this corresponds to the amplitudes linear in x and in Za, while for Delbriick scattering 
the amplitudes in the Born approximation are proportional to (Za) 2 . In Fig. [2] we plot the 
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ratio da /dan in units of % 2 /(Za) 2 . It can be seen that this ratio can be bigger than unity 
for the realistic parameters of the laser field. For example, for an optical (uq = 1 eV) laser 
with intensity of order of 10 23 W/cm 2 [ijj], for an initial photon energy iO\ = 500 MeV and 
for Z = 47 (silver), we obtain x 2 /(Za) 2 « 1.8. 




0.0 0.1 0.2 0.3 0.4 0.5 



P 

FIG. 2: The ratio da /dan in units x 2 /(Za) 2 obtained in the Born approximation. 



From Fig. [2] one sees that the region where the da mainly dominates over dan is that of 
small p. 

It is interesting to study the effect of higher order terms in x on the small- p differential 
cross section. This effect is determined by the Y-dependence of the function A(x) — 2|.Fi | 2 + 
I J~2 1 2 + I 1 2 (see Eq. (|35|) ). This dependence is shown in Fig. [31 The function A(x) has 
the following behaviour at small and large x : 

139 7904 2 

A(x) = 1 x at v < 1 , 

KAJ 4050 42525 A A 

, 3 1/3 5r 6 (2/3) s/ o 
^ = 4W378PW5) X ' atX>>1 ' (40) 

Therefore, the high order terms in x suppress, in comparison with the Born result, the cross 
section at large x an d increase it at small x-, making a substantial modification already at 
X<1- 



15 




FIG. 3: The ratio A(x)/A(0) , where A(0) = 139/4050. 



If the initial photon is not polarized, partial polarization of the final photon in the scat 
tering plane is produced. The corresponding Stokes parameters are [171] 

2E„Re(M n , ++ M; + _) 



^ 3 : E„(|M n , ++ | 2 + |M„, + _P) 
For p < 1 we find that S3 is p-independent and it is 



Si — S2 







S3 



2|J r i| 2 + LF 2 | 2 + l^sl 2 ' 



(41) 



(42) 



unlike the case of pure Delbruck scattering where S3 decreases at small p. The dependence 
of ^3 on x & t small p is shown in Fig. HI It can be seen that the polarization degree is not 
small and the x-dependence is smooth. 

It is also interesting to compare the cross section obtained here with that of photon 

nn 

splitting in a pure Coulomb field, see Refs. [5|, |6fl. The ratio of these two cross sections can 
be estimated as x 2 / a - Therefore, for realistic parameters of laser fields and photon sources 
this ratio can be substantially larger than unity. 

In conclusion, we have derived the Delbruck scattering amplitudes in a Coulomb field in 
the presence of a laser field. The amplitudes have been calculated in the Born approximation 
with respect to the Coulomb field and exactly in the parameters of the laser field having 
arbitrary strength, spectral content and polarization. The most interesting case of a high 
energy initial photon is investigated in detail for a monochromatic circularly polarized laser 
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0.54- 
0.52- 




0.0 0.5 1.0 1.5 2.0 

X 

FIG. 4: The Stokes parameter S3 at small p as a function of %. 



field. It is shown that the angular distribution of the final photon in the process under 
discussion substantially differs from that for Delbriick scattering in a pure Coulomb field. 
Moreover, the value of the cross section of the process at hand may exceed that of Delbriick 
scattering in a pure Coulomb field at realistic laser parameters and this essentially simplifies 
the possibility of the experimental observation of the phenomenon. Moreover, the effect of 
high order terms in x is quantitatively important already at relatively small Xi making the 
experimental investigation of these effects very attractive. 
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APPENDIX: COEFFICIENTS FOR THE HELICITY AMPLITUDES 



Here we present the coefficients -Bj> 1(T2 used in Eq. (J27J). We introduce the abbreviations 



Sj , (ri-r 2 )/ 2 . (ri+r 3 )/ 3 , (r 2 -ri)/ 2 (ri - r 3 )/ 3 

C(i — r : • (!■' - 



2r!r 2 



2r!r 3 



, ij2 + n)f 2 (n + r 3 )/ 3 7 , ? , j 

"3 = — ^ , n, = di + d 2 + d 3 



2t x t 2 2t x t 3 
(r 2 + 7i)/ 2 , (ri - r 3 )/ 3 



2riT 2 



2r!T 3 



+ 



2r!r 3 



2rir 2 

h h h 

g 1 = — - I - ih , g 2 = — - 1 - 2z/ 2 + iv 2 /i , # 3 = — - 1 - 2z/ 3 + ij/ 3 /i . (A.l) 



775 



r/S 



rjS 



Then, the nonzero coefficients for the contribution M^ ++ are 



B 2 ,++ = -4p £ 2 x x (g 2 g 3 + v 2 v 3 d\) , 

= -4£ 3 [g{g 2 g 3 + v 2 v 3 g\d\ + [y 3 g 2 + z/ 2 # 3 )|/i| 2 - %v 2 v 3 h*{h 2 + d 2 

2 [(8x1X3-1)5-2 . 4xix 3 # 3 . .^2, . , 
+p 4 1 h 1 — h + Aixidi 

L ^3 ^2 ^3 

4i£ £ 
+ ^Q?2 + gs) + -Q 2 {-4(t! + r 3 )(# 2 - tu 2 h)x 1 
a v 3 

— iiu 2 xi[2Tih + v 3 {2di + d 2 + d 3 )] + # 2 - iu 2 h} , 



5-1, 



4pf 



—9*92 + — gig* - x x g 2 gl + - g 2 d* 3 ) + h*di - v 2 hd* 3 

v 3 v 2 



^3 ^2 



p Aip 

xi(8x 2 x 3 - 1) - x 3 + — l xi 

v 2 v 3 S V ^3 ^2 



2x 2 x 3 



-i\g* x -ih*){g 2 -iu 2 h)+ p 



2X1X3 
^2 



r 2\ , 4z, 



1 Q 2 [ , 

^3 P I 

-4?V 2 £ 2 (^ - z/i*)[2n/i + z/ 3 (2di + d 2 + d 3 



-4^% 2 - iv 2 h)[2nh* + i/ 3 (dj + d*)] - p 2 l Xl u 3 + —)+ + 73) 

4£ 3 [-5i*525 3 + M^g*i - 9l)\h\ 2 + v 3 g 2 df - iv 2 v 3 h(h* 2 + d* 2 )] 

P 2 t ( A * 8x2X3-1 4z£ 

4xix 3 ^ 3 ^ + 4zx 3 d 3 -%—)+ — (g 1 + g 3 ) 

V 3 V 3 I D 



4 

^2 

+^g 2 
^3 

4p^3 

^2 



4x 3 (ri + T 3 ) 

^2 



1 — ^3 



( 5 *_^) + l!^[2 Tl / i * + i , 3 (d* + d*)] 

^2 



(A.2) 
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The nonzero coefficients for the contribution Ml 3 |_ are 



Bi, + - 



= p x 2 



ft 

3 8X1X3 — 1 



[(4X1X3 - 1)5-3 + ^3/1] , 



-4p^ 2 
Aip 



93 



V2 



X 3 

S \ v 2 



x\g* 2 ) - %3 - iv 3 h)d£ + v 3 ( ^ + x x v 2 ) \h\ 2 



At , \ XlX3 , , X3 1 

4(ri + r 3 ) h 1^1 H 1 

v 2 V2 



B-i,+- = 



P 2 i 



8x2X3 — 1 
^2^3 



-0i 



8x1X2 



^3 



1 # 4xix 3 - 1 „ 
-92 7. 93 



V2 



2i(l-2x 3 ) ^ | 2i{\ -4x 2 ) 



-2i(l - 2xiW* + 

^2 ^3 
-4£ 3 [-<?i*<? 3 <? 2 * + ^2<7i* - 9* 2 )\h\ 2 + V293df ~ iu 2 U 3 h(h* 2 + df)] 



4i£ 



-^(^! + ^) + -Q 2 -4(n + r 3 ) 

+ (1 + v 3 )(gl - ih*) - 4ix 1 z/ 2 [2r 1 /i* + v 3 (2d\ + d* + d* 3 )} 
Aix 3 



- Xi(# 2 + W2h*) 



-g 2 — iu 2 h* 



V2 



-finh* + v 3 (dl + d* 2 )\ 



+2i[2n(d* 2 + d* 3 ) + (2r 2 - v 2 ){d\ + d* 3 ) + 2r 3 K + d* 2 ) 



B_ 2 , + - = Apt 



■ i9*29t - 9*1 [ —9*2 + —9*3 ) + Wldl + ig* 2 d* 3 + ig* 3 d* 2 

v 3 v 2 



+x 1 is 2 is 3 d? + ^^d* 2 + —d? + v 2 d\d* 2 + v 3 d\d* 3 - d* 2 d* 3 



V3 



V2 



4f2 

+— Q 2 {(9*i ~ ih*)[(n + r 3 )(g* 2 + iis 2 h*) - 2iu 2 nh* - iv 2 v 3 {2d\ + d* 2 + d* 3 )} 

P v 3 

+i(g* 2 + ii/ 2 /i*)[2n/i* + v 3 (d\ + d* 2 )} + 2v 2 {d\ + d;)[n(d* 2 + d* 3 ) + T 3 {d\ + d* 2 )}} , 



3,H 



4£ {.9*i9*29l + v 2 v 3 g\dX - v 3 g* 2 d* 3 - u 2 g 3 d* 2 + 2iu 2 u 3 d* 1 d* 2 d* 3 ) 



(A.3) 



Finally, the coefficients for the amplitudes M^l_ + and M^!__ can be obtained with the help 
of the substitutions 



B jt -+ = B*_ jt+ _(S -> -S) , Bj-- = B*_ j++ (S -> -£) , 



(A.4) 
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where the replacement S — > —S means that after complex conjugation, it is also necessary 
to change the sign of the terms containing S in Eqs. (IA.2j) and flA.31) . 
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